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Estimators  for  model- based  passive 
localisation 

B  J  .Sullivan,  W.  Volkm&nn  and  S.  Bongi 


Executive  Summary!  1  stive  localiaation  techniques  have  not  been  re¬ 
ally  effective  for  many  yea  and  the  aim  of  thia  S  AC  LA  NIC  F  N  work  is 
an  attempt  to  improve  the  situation  Preliminary  results  or  thi  work  have 
shown,  for  certain  conditions,  a  radical  improvement  over  the  older  basic 
methods  bat  ranch  work  still  remains  to  bo  done  to  show  that  the  method 
remains  robust  in  a  variety  of  practical  et-aea  conditions. 

The  present  'standard'  techniques  of  bearinga-only  and  wavefront  curva¬ 
ture  make  so  use  of  available  environmental  information  about  the  medium. 
Measurement  of  the  multipath  arrivals  (rs-iging  on  the  vertical)  uses  some 
of  the  medium  information.  However,  the  newer  model-bused  (matched- 
held)  passive  ranging,  he  described  in  this  memorandum,  provides  better 
It  ablation  estimates  by  using  all  the  available  environmental  information. 

f  In  this  technique,  a  model  hi  selected  that  describes  the  propagation  channel 

to  a  reasonable  degree  and  thi*  is  used  to  make  a  prediction  of  the  eld 
received  at  an  array  This  in-  tel  «»W  <s  ja;s  compared  to  the  mess,  red 
nuSd  wu  a  see  of  source  coordinates  that  provides  the  best  match  is  taken 

t  as  an  estimate  of  the  source  position.  This  can  be  done  in  two  ways:  the 

problem  cun  be  inverted  and  one  simply  solves  for  the  source  coordinates  or 
a  search  can  be  made  over  a  prescribed  set  of  possible  source  coordinates. 
Bor  each  set  of  source  coordinates,  an  estimator  is  computed  and  used  to 
select  the  best  estimate  of  the  coordinate*. 

This  memorandum  is  a  ststi  repot  t  which  studies  the  effect  of  four  different 
types  of  estimator  using  both  synthetic  and  real  at-sea  data.  A  comparison  is 
made  between  Bucket’s  Estimator  and  thr>  a  types  of  least-squares  estimator. 
It  is  shown  that  Bttckor’u  Estimator  is,  in  general,  inferior  to  the  other  three 
but  it  is  recommended  that  further  work  with  at-sea  data  it  required  before 
one  eat-  comment  statistically  on  the  optimum  estimator 
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Estimator!  for  model- baaed  passive 
localisation 


E.J.  Sullivan,  W.  Volkraann  and  S.  Bongi 

Abstract:  A  comparative  study  is  made  of  the  performance  of  four  differ¬ 
ent  estimators  as  used  in  the  matched-Aeid  technique  of  passive  localisation. 
The  study  is  based  on  both  real  and  synthesised  data.  In  the  synthesised 
data  case,  a  comparison  is  made  of  the  performance  of  the  estimators  for 
various  signal-to-noise  ratios.  The  four  estimators  studied  are  Bucher’s  Es¬ 
timators,  which  can  be  thought  of  as  a  spatial  matched  Alter,  and  three 
likelihood-type  estimators.  The  results  indicate  that  the  matched-Aeld  type 
estimator  has  a  slightly  better  signal-to-noiae  performance  than  the  others, 
but  rather  poor  eideiobe  behaviour,  whereas  for  the  likelihood-type  estima¬ 
tor  the  sidelobe  behaviour  is  quite  good.  N 


Keywords:  inverse  problentSjO  matched-Aeld  processing  o  model-based 
signal  processing  •’^passive  localisation  w1  passive  ranging 
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1.  Introduction 

Model-based  passive  ranging,  sometimes  referred  to  as  matched-field  processing,  is 
basically  an  inverse  problem.  A  propagation  model  that  describes  the  propagation 
channel  in  some  acceptable  sense  is  selected,  and  is  used  to  make  a  prediction  of 
the  field  received  on  an  array  (usually,  but  not  necessarily  vertical).  This  model 
field  is  then  compared  to  the  measured  field  and  the  set  of  source  coordinates  that 
provide  the  best  match  to  the  measured  field  is  then  taken  as  the  estimate  of  the 
source  coordinates.  This  procedure  can  be  carried  out  in  two  different  ways.  The 
problem  can  be  directly  inverted,  solving  for  the  source  coordinates,  or  a  search  can 
be  made  over  a  prescribed  set  of  source  coordinates  in  some  manner.  For  each  set  of 
source  coordinates,  an  estimator  is  computed  and  used  to  select  the  best  estimate 
of  the  coordinates.  In  the  case  of  range-depth  estimation,  this  estimator  can  be 
plotted  on  a  range-depth  map  such  that  the  estimates  can  be  directly  taken  as  the 
coordinates  of  some  extremum.  It  is  this  technique  that  we  are  concerned  with  here 
In  particular,  we  make  a  comparative  study  of  four  different  estimators  based  on  both 
real  and  synthetic  data.  The  estimators  are  Bucker’s  Estimator  and  three  types  of 
least-squares  estimator.  A  thorough  discussion  of  the  subject  of  model-baaed  passive 
localisation  can  be  found  in  [1], 
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2.  Theory 


For  the  layered  waveguide  model  with  source  on  tne  z  (or  vertical)  axis,  the  pressure 
field  is  symmetric  about  z  and  is  therefore  governed  by  the  cylindrical  wave  equation 
which  is  given  by 


|^P(r,r1t)+  i|:p(r1r1t)+  z,t)=  ^X-p(r,z,t). 


1 


(1) 


Since  in  this  study  a  shallow-water  range-independent  scenario  is  chosen,  Eq.  (1)  can 
be  solved  by  separation  of  variables,  since  the  range  independence  allows  a  factored 
form  for  the  solutions.  Assuming  far-field  conditions  and  a  harmonic  source  (see  [2] 
for  details)  the  normal-mode  model  of  Eq.  (2)  obtains: 

p(r- z)  =  X]  7E=ei*,r'  (2) 

m=l  V*"»P 


Here  p  is  the  acoustic  pressure,  <f>m(z)  is  the  mth  modal  function  evaluated  at  r,  zo 
is  the  source  depth,  r  is  the  horizontsd  range,  km  is  the  horizontal  wave  number  for 
the  mth  mode,  am  is  the  loss  factor  for  the  mth  mode,  and  M  is  the  number  of 
modes. 


Assuming  a  vertical  receiving  array,  Eq.  (2)  can  be  conveniently  written  in  matrix 
form  as 


Pti  —  Afn,mX'ni 


(3) 


where  Pn  is  the  pressure  at  the  nth  hydrophone  of  the  vertical  receiving  array  with 
N  hydrophones,  and 


"  VK  ’ 


Xm  — 


^(zo)*-*0"-^ 


The  estimation  algorithm  proceeds  as  follows: 

(1)  Given  the  sound- velocity  profile  (SVP),  the  ocean  depth  and  the  bottom 
boundary  conditions,  compute  <t>m(zn),  ^m(i0),  km  and  om  for  all  desired 
values  of  z® 

(2)  Compute  the  model  pressure  vector  {F^(ro>  r)}  for  a  given  set  of  values  of 
r  and  x«,  from  Eq.  (3). 
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(3)  Compare  {Pi1}  to  the  data  vector  [Pj?}  over  the  prescribed  ranges  of  zg 
and  r. 

It  is  the  third  step  that  we  are  concerned  with  here.  Historically  the  estimator  used 
for  step  3  for  the  forward  modelling  procedure  was  based  on  the  inner  product  of 
{ and  {P®}  (3).  This  we  will  refer  to  as  Bucher’s  Estimator.  The  form  we  use 
is  written  as 

i>bu  =  OV(I-pT  |j,jT),  «) 

m 

where  |P'D|J  and  |PMjJ  are  the  squared  magnitudes  of  the  data  and  model  vectors, 
respectively.  We  note  that  if  one  computes  the  expected  value  of  £>bu,  Bq.  (4) 
becomes 

l  m 

=  (6) 

1  m 

where  flm,i  is  the  covariance  matrix  of  the  data.  In  Eq.  (5)  the  normalization  has 
been  ignored  for  simplicity;  however  in  this  study  we  will  use  the  form  expressed  by 
Eq.  (4).  One  then  seeks  the  values  of  r  and  io  that  maximize  Dbv 

The  second  estimator  that  we  consider  is  the  least  squares  fit  of  PM  to  PD ,  i.e.  we 
seek  the  minimum  of 

t>LS=  ZlP*-P£l’-  (6) 

m=  1 

The  data  and  model  vectors  are  normalised  to  a  reference  hydrophone. 

The  third  estimator  is  the  incoherent  form  of  Di s.  That  is,  we  fit  the  hydrophone 
powers  in  a  least-squares  sense.  Thus  we  seek  the  minimum  of 

0PW=£(l^l’-lW-  m 

m=l 


The  fourth  estimator  is  based  on  an  inversion  of  Eq.  (3),  and  since  this  equation 
is  linear  in  \m,  it  can  easily  be  solved  by  the  method  of  least  squares,  under  the 
assumption  that  N  >  hi.  This  results  in 

X,  =  [(M'M)-'Af'],.nPn,  (8) 
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where  M'  is  the  transpose  of  M .  The  expression  (M1  M)~l if1  is  sometimes  referred 
to  as  the  pseudo-inverse  of  M.  Naturally,  this  reduces  to  M~l  when  N  =  M. 
Equation  (8)  actually  constitutes  a  modal  filter,  thus  allowing  modal  selection  in 
the  solution.  This  turns  out  to  be  important  in  our  case.  The  estimator  we  use  for 
this  approach  is  based  on  the  least-squares  fit  of  xi  to  the  prediction  of  xi  by  the 
model.  That  is,  we  seek  the  minimum  of 

0ml  =  5>,"-X,T,  <•> 

*=1 

where  ML  refers  to  maximum-likelihood  in  order  to  differentiate  this  estimator  from 
fits,  and  K  refers  to  the  number  of  modes  used.  Here,  os  in  the  case  of  Zips,  the 
model  and  data  vectors  are  normalised  to  a  given  complex  amplitude. 
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3.  Data  characteristics 

As  mentioned  in  the  introduction  both  real  and  synthetic  data  were  used  in  this 
study.  The  synthetic  data  were  generated  by  the  SNAP  model  [4],  An  ocean  depth 
of  103  m  was  selected  with  a  winter  Mediterranean  sound  velocity  profile  (SVP).  A 
point  acoustic  source  of  190  Hz  was  taken  to  be  at  a  depth  of  50  m  and  a  range  of 
7.4  km  from  the  receiving  array.  This  situation  supported  9  modes.  The  receiver 
was  a  vertical  array  of  32  point  hydrophones  with  a  spacing  of  2  m.  The  topmost 
hydrophone  was  located  at  a  depth  of  30  m.  The  bottom  was  composed  of  2.5  m 
of  sand  over  a  rock  sub-bottom.  These  complex  synthetic  data  were  then  modified 
by  the  addition  of  a  complex  noise  term  that  was  taken  from  a  random  number 
generator  with  gaussian  statistics.  The  signal- to- noise  ratio  at  each  hydrophone  was 
taken  to  be  equal. 

The  Teal  data  were  taken  at  a  depth  of  103  m  in  the  region  north  of  the  island  of 
Elba  in  the  Mediterranean.  The  SVP  and  bottom  conditions  were  those  used  for  the 
generation  of  the  synthetic  data.  The  source  and  receiver  coordinates,  as  well  as  the 
frequency  were  also  the  same  as  for  the  synthetic  data.  The  data  were  preprocessed 
by  performing  a  256-point  FFT  on  a  2  s  data  record.  The  frequency  line  at  190  Hz 
was  then  taken  as  a  ‘snapshot’  of  the  complex  amplitude.  The  S/N  ratio  of  these 
data  was  30-40  dB. 

Due  to  environment  ail  conditions,  mainly  currents,  the  vertical  array  was  not  always 
vertical.  Hence,  it  must  be  assumed  that  in  the  case  of  the  real  data,  there  are  errors 
in  the  assumed  positions  of  the  hydrophones  in  the  horizontal  direction. 
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4.  Results 

The  four  estimators  were  studied  with  regard  to  their  performance  under  various 
S/N  ratios  in  the  case  of  the  synthetic  data,  and  also  their  performance  with  real 
data.  Figures  1  through  3  depict  the  performance  of  Vbv  for  S/N  ratios  of  10,  20 
and  30  dB  respectively.  Figure  4  gives  the  results  of  D bu  for  real  data.  Figures 
5  through  8  depict  the  same  series  of  cases,  i.e.  10,  20,  30  dB  and  real  data,  for 
2?ls-  Continuing  in  this  manner,  Figs.  9  through  12  show  the  same  four  cases, 
respectively,  for  Z?pw  and  finally,  Figs.  13  through  16  present  these  same  four  cases 
for  £>mL‘  As  mentioned  in  Sect.  2,  we  seek  the  maximum  of  £>bv  and  the  minimum 
of  the  other  three  estimators.  However,  in  these  illustrations,  the  likelihood  surface 
has  been  inverted  for  Du,  Dm  and  Dp  w  so  that  the  extremum  that  we  seek  is  the 
maximum  in  all  four  cases.  This  was  done  by  dividing  all  values  of  the  likelihood 
surface  into  its  minimum  value  for  these  cases. 

There  are  severed  conclusions  that  can  be  drawn  from  these  figures.  First,  we  see 
that  Dgu  suffers  from  a  severe  sidelobe  problem,  while  at  the  same  time  has  the 
best  S/N  performance.  Although  difficult  to  see  due  to  the  clutter  caused  by  the 
sidelobes,  the  maximum  of  £>bu  occurs  at  the  correct  position  of  r  =  7.4  im  and 
d  ~  50  m.  In  fact,  this  behaviour  holds  down  to  a  S/N  ratio  of  0  dB  (which  is  not 
shown).  We  note,  however,  that  Dbv  failed  to  indicate  the  solution  in  the  case  of 
the  real  data. 

Considering  next  the  behaviour  of  Dis  >  we  see  that  there  is  an  immediate  improve¬ 
ment  in  sidelobe  performance,  but  the  S/N  performance  is  not  quite  as  good  since 
Dis  failed  to  provide  a  solution  at  S/N  =  0  dB  (which  case  is  not  shown).  We  also 
see  that  ZJLS  does  not  provide  a  solution  for  the  real  data  case. 

Continuing  on  to  Upwi  we  find  the  S/N  performance  slightly  degraded  with  respect 
to  that  of  Du ,  but  there  is  a  solution,  although  badly  aliased,  in  the  case  of  the  real 
data.  Although  not  obvious  in  the  illustration,  the  peak  at  the  correct  coordinates 
of  r  =  7.4  km  and  d  =  50  m  is  essentially  equal  to  those  at  13.6  km,  d  =  60  m  and 
r  =  5.0  km,  d  —  15  m.  Thus,  the  solution  is  not  unique.  It  should  be  noted  that  the 
solution  seems  to  be  unbiased. 

Proceeding  to  the  last  case,  that  of  X?ml,  we  find  that  the  S/N  performance  is  even 
nuTac  that,  that  of  D pw  •  Xlovfewi,  again  there  is  a  solution  witn  the  real  data.  It  is 
biased  in  depth  and,  as  can  be  clearly  seen,  there  are  two  larger  erroneous  solutions. 
This  solution  was  obtained  by  using  only  the  first  7  modes  in  the  solution.  If  more 
or  less  than  7  modes  were  used,  no  solution  was  obtained. 
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5.  Discussion 

Generally  speaking  the  performance  of  Dbu  seems  to  be  inferior  to  the  other  three 
estimators.  Although  its  S/N  performance  is  better  than  that  of  any  of  the  others, 
the  sidelobe  behaviour  is  so  poor,  that  without  a  perfect  match  to  the  model  one 
would  expect  this  advantage  to  disappear  rapidly,  since  a  slight  elevation  of  any  of 
the  aliased  solutions  would  be  disastrous.  The  sidelobe  behaviour  could  be  improved 
in  two  situations.  First  more  modes  tend  to  lower  the  sidelobes;  therefore  a  scenario 
with  higher  frequency  and/or  deeper  ocean  depth  could  improve  things.  Secondly, 
the  array  data  can  be  ‘beamformed’,  i.e.  can  be  passed  through  some  preprocessing 
to  provide  more  spatial  selectivity.  This  is  done  in  [4]  where  the  so-called  maximum- 
likelihood  beamformer  is  used  to  preprocess  the  data. 

Since  the  hydrophone  positions  were  not  well  known  due  to  the  tilting  of  the  vertical 
array,  one  suspects  that  the  generally  poor  performance  in  the  case  of  the  real  data 
is  due  to  this  ‘mismatch’  between  the  model  and  the  true  situation.  Thus  one 
could  speculate  that  the  phase  errors  produced  by  this  mismatch  are  ignored  by 
Dpw,  which  is  essentially  an  incoherent  version  of  Dj, s.  Thus  at  the  expense  of 
some  S/N  performance,  robustness  to  mismatch  is  achieved  by  Dpw .  This  issue  of 
mismatch  enters  into  the  case  of  Dml  also.  Here,  as  pointed  out  above,  a  solution 
was  obtainable  only  in  the  case  of  7  modes.  This  suggests  that  the  mismatch  caused 
unacceptable  errors  in  modes  8  and  9  and  therefore  the  ability  to  filter  them  out  is 
crucial  to  a  coherent-type  estimator. 

A  general  conclusion  then,  is  that  array  tilt  is  a  serious  problem  and  the  level  of 
seriousness  can  depend  on  the  particular  estimator  that  one  uses.  One  way  that 
such  problems  can  be  avoided  is  to  eliminate  the  need  for  a  priori  knowledge  of 
the  model  parameters,  A  means  by  which  this  can  be  done  is  given  in  [5]  where 
it  is  shown  that  knowledge  of  the  horizontal  wavenumbers,  which  can  be  estimated 
with  either  a  long  towed  array  or  a  synthetic  aperture  towed  array,  is  sufficient  to 
estimate  range.  This  could  have  important  ramifications  for  these  methods,  since 
a  recent  study  [6]  has  shown  that  mismatch  in  the  environmental  parameters  can 
have  serious  deleterious  effects. 

Finally,  it  is  worth  mentioning  that  although  the  general  S/N  performance  does  not 
seem  too  promising,  the  synthetic  data  were  not  averaged.  That  is,  the  noise  was 
directly  added  as  a  single  realization  of  a  gaussian  process.  This  corresponds  to  a 
single  ‘snapshot1  of  real  data  with  no  time  averaging.  Thus  the  S/N  values  should 
be  considered  as  a  worst- case  situation  that  would  improve  with  averaging. 
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Fig.  1.  Contour  plot  of  Bucker’s  Estimator  for  the  case  of  synthetic  data  with  a  signal- to- 
noise  ratio  of  10  dB. 
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Fig.  3.  Contour  plot  of  Bucket's  Estimator  for  the  case  of  synthetic  data  with  a  signal-to- 
noise  ratio  of  20  dB. 
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Fig .  5.  Contour  plot  of  the  mean-squared  error  of  tke  pressure  field  for  the  case  of  10  dB 
signal- to- noise  ratio.  The  data  are  inverted  sack  that  tke  maximum  of  the  plot  corresponds 
to  the  minimum  mean- squared  error. 
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Hr.  «.  Contour  plot  of  the  meu-aqaared  error  of  tke  preeeute  field  for  the  cue  of  30  dB 
elgnei-to-noise  r»tk>.  The  date  are  inverted  inch  that  the  maximum  of  the  plot  correspond* 
to  the  minimum  uiean-<qoered  error. 


-  03  - 


S  ACL  ANTCKN  SU-11 1 


□ 


> 

O.M 

O.M 

O.M 

O.M 

0.10 


1.00 


-  O.M 

-  OM 

-  O.M 

-  O.M 
0.30 


Level 


Fig.  7.  Contour  plot  of  the  mean-equated  error  of  the  pressure  field  for  the  case  of  30  dB 
sign al-to- noise  ratio.  The  data  are  inverted  such  that  the  maximum  of  the  plot  corresponds 
to  the  minimum  mean-squared  error. 
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Fi  t  «  Contour  plot  of  the  mean-squared  eiroi  of  the  pressure  field  for  the  ewe  of  teal  date. 
The  data  ate  inverted  such  that  the  maximum  of  the  plot  corresponds  to  the  minimum 
mean-squared  error. 
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Fig.  9.  Contour  plot  of  the  mean  squared  error  of  the  hydrophone  power  for  the  case 
of  synthetic  data  with  signal- to-noise  ratio  of  10  dB.  The  data  are  inverted  such  that  the 
maximum  of  the  plot  corresponds  to  the  minimum  mean-squared  error. 
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fig.  10.  Contour  plol  of  the  mean-squared  error  of  the  hydrophone  power  for  the  cue 
of  eynthetk  dit>  with  signal-to-noise  ratio  of  30  dB.  The  date  ue  inverted  such  that  the 
maximum  of  the  plot  correspond*  to  the  minimum  mean-squired  error. 
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Fig.  11.  Contoui  plot  of  the  mean-squared  eiror  of  the  hydrophone  power  for  the  cue 
of  synthetic  date  with  signal-to-noise  ratio  of  20  dB.  The  data  ate  inverted  such  that  the 
maximum  of  the  plot  correspond*  to  the  minimum  mean-equated  ettot. 


Fig ■  12.  Content  plot  of  the  mean-squared  error  of  the  hydrophone  power  for  the  case  of 
teal  data.  The  data  ate  inverted  such  that  the  maximum  of  the  plot  corresponds  to  the 
minimum  mean-squared  ettor. 
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Fig.  13.  Contour  plot  of  the  mean-squared  error  of  the  modal  amplitudes  for  the  case 
of  synthetic  data  vith  signal- to- noise  ratio  of  10  dB.  The  data  are  inverted  such  that  the 
maximum  of  the  plot  corresponds  to  the  minimum  mean-squared  error. 


Fig.  14.  Contour  plot  of  the  mean-squared  error  of  the  modal  amplitudes  for  the  case 
of  synthetic  data  with  si(nal-to-noise  ratio  of  20  dB.  The  data  are  inverted  such  that  the 
maximum  of  the  plot  corresponds  to  the  minimum  mean-squared  error. 
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Fig-  IS.  Contour  plot  of  the  mean-squared  error  of  the  model  amplitude*  for  the  cnse 
of  synthetic  date  with  sijnal-to-noise  ratio  of  30  dB.  The  data  are  inverted  etich  that  the 
maximum  of  the  plot  correspond*  to  the  minimum  mean-*quared  error. 
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Fig.  16.  Contour  plot  of  the  mean-squared  error  of  the  modal  amplitude*  for  the  case  of 
real  data.  The  data  are  inverted  such  that  the  maximum  of  the  plot  corresponds  to  the 
minimum  mean-squared  error. 
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